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Abstract 

It is known that strong electric fields produce electron and positron pairs from the 
vacuum, and due to the back-reaction these pairs oscillate back and forth coherently 
with the alternating electric fields in time. We study this phenomenon in spatially 
inhomogeneous and bound electric fields by integrating the equations of energy- 
momentum and particle-number conservations and Maxwell equations. The space 
and time evolutions of the pair-induced electric field, electric charge- and current- 
densities are calculated. The results show that non- vanishing electric charge-density 
and the propagation of pair-induced electric fields, differently from the case of homo- 
geneous and unbound electric fields. The space and time variations of pair-induced 
electric charges and currents emit an electromagnetic radiation. We obtain the nar- 
row spectrum and intensity of this radiation, whose peak cu pea k locates in the region 
around 4 KeV for electric field strength ~ E c . We discuss their relevances to both 
the laboratory experiments for electron and positron pair-productions and the as- 
trophysical observations of compact stars with an electromagnetic structure. 
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Introduction. As reviewed in the recent report pQ, since the pioneer works by 
Sauter j2], Heisenberg and Euler [3] in 1930's, then by Schwinger jl] in 1950's, 
it has been well known that positron-electron pairs are produced from the 
vacuum in external electric fields. In a constant electric field Eq in dependent 
of space and time, the pair-creation rate per unit volume is given by [3f4"] . 

dN m\ fE \ 2 f E c \ 
S = = TZs Ur exp -tt— , (1) 



dVdt A7r 3 \Ej * V E, 



-0 



where the critical field E c = m 2 e c'/{eh), the Plank's constant H, the speed of 
light c, the electron mass m e , the absolute value of electron charge e and the 
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fine structure constant a = e 2 /he (in this article we use the natural units h = 
c = 1, unless otherwise specified). The pair-production rate ([1]) is significantly 
large for strong electric fields E > E c ~ 1.3 • 10 16 V/cm. The critical field 
will probably be reached by recent advanced laser technologies in laboratory 
experiments [5|6|7] , X-ray free electron laser (XFEL) facilities [8] , optical high- 
intensity laser facilities such as Vulcan or ELI [9J, and SLAC E144 using 
nonlinear Compton scattering [10] . On the other hand, strong overcritical 
electric fields (E > 10 E c ) can be created in astrophysical environments, for 
instance, quark stars [TT1fr2] and neutron stars [T3]-[rB]. 

The back-reaction and screening effects of electron and positron pairs on exter- 
nal electric fields lead to the phenomenon of plasma oscillations: electrons and 
positrons moving back and forth coherently with alternating electric fields. 
This means that external electric fields are not eliminated within the Comp- 
ton time h/m e c 2 of pair-production process, rather oscillate collectively with 
the motion of pairs in a much longer timescale. 

In a constant electric field E (JT]), the phenomenon of plasma oscillations is 
studied in the two frameworks PQ: (1) the semi-classical QED with quantized 
Dirac field and classical electric field [TTIfTS] ; (2) the kinetic description using 
the Boltzmann-Vlasov and Maxwell equations |19}20B21| 22|23j. In the second 
framework, the Boltzmann-Vlasov equation is used to obtain the equations 
for the continuity and energy- momentum conservations [20J. 

Ref. |23J shows the evidence of plasma oscillation in under-critical field {E < 
E c ) and the relation between the kinetic energy and numbers of oscillating 
pairs in a given electric field strength Eq. Taking into account the creation and 
annihilation process e + +e~ <^ 7+7, it is shown [22] that the plasma oscillation 
in an overcritical field is led to a plasma of photons, electrons and positions 
with the equipartition of their number- and energy-densities. The phenomenon 
of plasma oscillations is studied in connection with pair creation in heavy ions 
collisions [19]- [21], the laser field [21]- [27], and gravitational collapse [28]. It 
is worthwhile to emphasize that the plasma oscillation occurs not only at 
overcritical field-strengths E > E c (see for instance Refs. [T7f22] ). but also 
undercritical field-strengths E < E c (see Ref. [23]), and plasma oscillation 
frequency is related to field-strength E , while the number of oscillating pairs 
depends on the pair-production rate (Q. More details can be found in the 
recent review article pp. 

The realistic ultra-strong electric fields are not only vary with space and time, 
but also confined in a finite region. In this letter, studying the plasma os- 
cillations in spatially inhomogeneous electric field, we present the evidence 
of electric fields propagation, leading to electromagnetic radiation with a pe- 
culiar narrow spectrum in the KeV-region, which should be distinctive and 
experimentally observable. 
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In the kinetic description for the plasma fluids of positrons (+) or electrons 
(— ), whose single-particle spectrum p\. = (p^ + ml) 1 / 2 , we define the number- 
densities n±(t,x) and "averaged" velocities v±(t,x) of the fluids: 

n±(t,x) = J ^^/±(t,P±,x), (2) 

V±(t ' x) "^/{2#(t) /±(t ' P± ' X) ' (3) 

where /±(t, p±,x) is the distribution function in the phase space. The four- 
velocities of the electron and positron fluids U± = 7±(1, v±), the Lorentz fac- 
tor 7± = (1 — |v±| 2 ) -1 / 2 , and the comoving number-densities n± = n±( r y±)~ 1 , 
where we choose the laboratory frame where pairs are created at rest. The 
collision-less plasma fluid of electrons and positrons coupling to electromag- 
netic fields is governed by the continuity, energy-momentum conservation and 
Maxwell equations: 

(4) 



dT. 



/// 

± 



0r „ -F^ + J^J, (5) 

— = -M^Ld + ^pola + JL), (6) 

where S = dN/dVdt is the pair-production rate, J± = ±en±U± electric cur- 
rents and the energy-momentum tensors [30] 



rr = P±!T + (P± + e±)U£U£, (7) 

and the pressure p± and comoving energy- density e± is related by the equation 
of state, in general < p± < e±/3. In the laboratory frame, the fluid energy- 
density e± = T 00 and momentum-density p± = T l ° are given by 

e± = (e± + P±v 2 t )7|, p ± = (e± + P±)t±v±. (8) 

In Eqs. (ISfBl) F£ is the tensor of electromagnetic fields (E, B), the conducting 
four- current density 

^ ond = e(n + ^-n_^), d,J? ond = 0, (9) 
and polarized four-current density J£ ola = £ ± J± po \ a and ^± P oia = (Ppol 



± 

a' ° pola 



(2vr)3pO J 



^pda = K = I tt^St^ ( 10 ) 
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and S = J d 3 p±/[(2ii) 3 p±}S. Using "averaged" velocities ([3]) of the fluids, we 
approximately have 



± „ m el ±S t ± ^ m e7± \v ± \S 

J pola — ijji ■ Cj ' Fpola — ^ > 1111 



where the magnetic field B = 0. In Eq. (jSJ), J^ t = (/w, Jext) is an external 
electric current. 

Basic equations of motion. For simplicity to start with, we consider the 
electric field E ext created by a capacitor made of two parallel plates, one 
carries an external charge +Q and another —Q. The sizes of two parallel 
plats are L x and L y , which are much larger than their separation £ in the 
z-direction, i.e., L x ^> £ and L y ^> £. For \z\ ~ 0(£), the system has an 
approximate translation symmetry in the (x, y) plane. As results the electric 
field E ext (x, y,z) ~ E ext (z)z and B ext (:c, y, z) m 0, is approximately homo- 
geneous in the (x, y) plane and confined within the capacitor. In addition, 
<9E cxt /<9t Rj 0, namely, this electric field is assumed to be continuously sup- 
plied by an external source (+Q, —Q) or slowly varying. 

In order to do calculations we model this electric field as the one- dimensional 
Sauter electric field in the z-direction 



E ext (z) = E /cosh\z/£), a = eE £/m e c 2 = (£/ X C )(E / E c ), (12) 

where the Ac* is Compton wavelength, the external electric charge is given by 
dE ext (z)/dz = 47rpext and the external electric current vanishes J ext = for 
the field being static dE ext /dt = 0. In the electric field configuration (TT2"]) and 
B rj 0, the "averaged" velocities v± of electrons and positrons fluids are in 
the z-direction, 

L/£ = 7 ± (l,0,0,±u ± ), (13) 
and the total fluid current- and charge-densities (jBJ) «/ M = (p, J) are 



T , , m e ( 7+ + 7_)g 
J z = en + v + + en-V- H — , (14) 

p = e{n+ - n _ ) + ^±p^ll. (15) 

The system can be approximately treated as a 1 + 1 dimensional system in 
terms of space-time variables (z,t), and Eqs. (jMSJ) become for zero pressure 
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dn± dn±v± 

-w ± — = s ' (16) 

^ ± = en±v±E + m e 7±S', (17) 
dp± dp±v± 

-ST ± — ^ = en±E + m e 7 ± f ± 6, (18) 

dE , s 

^ = -4vrJ z , (19) 

_ = 47r(p + p ext ). (20) 

The total electric field E(z,t) in Eqs. (I16H19P is the superposition of two 
components: 

E(z,t) = E cxt (z) + E- md (z,t), (21) 

where the space- and time-dependent E ind (z,t) is the electric field created by 

electron and positron pairs. We call J z (z,t) f|T4|) . p(z,t) (fT5|) and E- md (z,t) 
pair-induced electric current, charge and field. 



As for the pair-production rate S in Eqs. (U6H19p . instead of the pair-production 
rate §Q for a constant field E , we adopt the following z-dependent formula 
for the pair-production rate in the Sauter field (TT21) . obtained by using the 
WKB-method to calculate the probability of quantum-mechanical tunneling 
E9], 



= m| E E(z) e _„ GmEc/E 
V ' ^E 2 G{0,£] V ; 

where G(0, £) and (5(0, £) are functions of the energy-level crossings £(z) and 
we approximately adopt E(z) « E Q /G(0,£) « E /G(0,£) in Eq. (122]) in order 
to do feasible numerical calculations. As shown by the Fig. 2 in Ref. [29] . the 
deviation of the pair-production rate ( 122]) due to this approximation is small. 
The formula ([22]) is derived for the static Sauter field (I12p . However, analo- 
gously to the discussions for the plasma oscillations in spatially homogeneous 
fields [22]- [21], it can be approximately used for a time- varying electric field 
E(z,t) (12T]) . provided the time-dependent component £^(,2, t), created by 
electron-positron pair-oscillations, varies much slowly compared with the rate 
of electron-positron pair-productions 0{m e c 2 /h). This can be justified by the 
inverse adiabaticity parameter f33]-|36j. 

n^e En , \ 

" = 17i >>1 - < 23) 

where u is the frequency of pair-oscillations. 
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Eqs. ( TT6|17||18p describe the motion of electron-positron plasma coupling to 
the electric field E and source S of pair-productions. The Maxwell equations 
(I19|20p describe the motion of the electric field ( ]2"T|) coupled to the current- 
and charge-densities (fT5]l . leading to the wave equation of the propagating 
electric field E [n d(z,t) [37J, 

d 2 E ind 1 d 2 E ind (?P,L!Mz\ (2A) 
dt 2 c 2 dz 2 K \dz c 2 dt / ' 1 j 

where we use dE ext /dz = 47rp ext and dE ext /dt = 0. This wave equation shows 
the propagating electric field E ind (z, t) in the region 1Z where the non-vanishing 
current J z and charge p are, and both the propagation and polarization of 
the electric field are in the z-direction. This implies a wave transportation of 
electromagnetic energies inside the region TZ. Since the current- and charge- 
densities (p, J z ) are functions of the field E(t,z) (ED) , the wave equation is 
highly nonlinear, the dispersion relation of the field is very complex and the 
velocity of field-propagation is not the speed of light. 

Numerical integrations. Given the parameters E = E c and i = 10 5 Ac of 
the Sauter field (|12|) as an initial electric field E ext , we numerically integrate 
Eqs. fU6tiTU]) in the spatial region TZ: —£/2 < z < £/2 and time interval T: 
< t < 3500T(7, where tq is the Compton time. The value T < 3500rc 
is chosen so that the adiabatic condition (|23|) is satisfied, and the spatial 
range TZ is determined by the capacity of computer for numerical calculations. 
The electric field strength Eq is chosen around the critical value E c , so that 
the semiclassical pair-production rate (122j) can be approximately used. Actu- 
ally, E ,£ and T are attributed to the characteristics of external ultra-strong 
electric fields E ext established by either experimental setups or astrophysical 
conditions. 

In Figs. [Hand El we respectively plot the time- and space-evolution of the total 
electric fields E(z,t) (I2T1) as functions of t and z at three different spatial 
points and times. As discussed in Figure captions, numerical results show 
the properties of the electric field wave E ind (z, t) propagating in the plasma of 
oscillating electron-positron pairs, as described by the wave equation (I21jl . This 
electric field wave propagates along the directions in which external electric 
field-strength decreases. The wave propagation is rather complex, depending 
on the space and time variations of the net charge density p(z, t) and current 
density j z (z, t), as shown in Figs. HH5J The net charge density p oscillates (see 
Figs. [3]and 0J proportionally to the field-gradient PU|) and at the center z = 
the charge density and field-gradient are zero independent of time evolution 
(see Fig.HJ). However, the total charge of pairs Q = J n d 3 xp must be zero at any 
time, as required by the neutrality. The electric current j z (z,t) alternating in 
space and time follows the space and time evolution of the electric field E(z, t) 
see Eq. (TH?j) . as shown in Figs. \5\ and El 
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We recall the discussions of the plasma oscillations in the case of spatially 
homogeneous electric field E without boundary [22 .23J . Due to the spatial 
homogeneity of electric fields and pair-production rate S ([I]), the number- 
densities n±(t,x) = n(t) (PD, "averaged" velocities \v±(t, x)| = v(t) (J3J) and 
energy-momenta e±(£, x) = e(t), |p±(t, x)| = p(t) (JBJ) are spatially homoge- 
neous so that the charge density ffT5]) p = identically vanishes and current 
CHD Jz = J z (t). All spatial derivative terms in Eqs. ffToTfTB]) and Eq. (PJU) 
vanish and Eq. (120]) becomes irrelevant. As results, the plasma oscillations 
described is the oscillations of electric fields and currents with respect time 
at each spatial point, and the electric field has no any spatial correlation and 
does not propagate. 

In contrary to the plasma oscillation in homogeneous fields, the presence of 
such field-propagation in inhomogeneous fields is due to: (i) non-vanishing 
field-gradient d z E (1201) and net charge- density p (JT5l) . as shown in Figs. [3] and HJ 
give the spatial correlations of the fields at neighboring points; (ii) the stronger 
field-strength, the larger field-oscillation frequency is, as shown in Fig. [TJ (iii) 
at the center z = the field-strength is largest and the field-oscillation is 
most rapid, and the field-oscillations at points \z\ > are slower and in retard 
phases, as shown in Fig. [2j The point (i) is essential, the charge density p 
oscillates (see Figs. [3] and HJ) proportionally to the field-gradient Eq. f[2"U|) and 
at the center z = the charge density and field-gradient are zero independent 
of time evolution (see Fig. HJ). Such field-propagation is reminiscent of the drift 
motion of particles driven by a field-gradient ( "ponderomotive" ) force, which 
is a cycle-averaged force on a charged particle in a spatially inhomogeneous 
oscillating electromagnetic field 



Radiation fields. As numerically shown in Fig. UM\ the propagation of the 
electric field wave E ind (z,t) inside the region 71 is rather complex, due to 
th high non-linearity of wave equation (JMJ). Nevertheless, the electromagnetic 
radiation fields E ra d and B rac j far away from the region 1Z are completely deter- 
mined and could be experimentally observable. At the space-time point (t, x) 
of an observer, the electromagnetic radiation fields E rad (;z, t) and B rad (z,i), 
emitted by the variations of electric charge density p(x', t') and current-density 
J(x',t') in the region 1Z (x' G 1Z) and time t' (t' G T), are given by 



E r ad(t,x)=-^rf 3 x'{— [ P (t',x')] 

~dJ{t',x' 



R 

ret + cR 



+ 



c 2 R 



dt' 



- ret 



B rad (t,x) =^rf 3 x'|[J( 



R 



dp(t',x') 
dt' 



(t',x')] ret x^ + 



aj(f,x') 

dt' 



ret 



R 

~(Fr 



(25) 
(26) 



where the subscript "ret" indicates t' = t — R/c, R = |x — x'|. In the radiation 
zone |x| 3> |x'| and R ~ |x|, where is far away from the plasma oscillation 
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Fig. 1. Electric fields E(z,t) are plotted as functions of t at three different points: 
z = (red), z = l/A (blue) and z = 1/2 (black). Analogously to the plasma 
oscillation in homogeneous fields, the stronger initial field-strength, the larger field- 
-oscillation frequency is, i.e., ui(z = 0) > uj(z = £/A) > u)(z = £/2), where u)(z) is 
the field oscillating frequency at the spatial point z. 



Fig. 2. Electric fields E(z,t) are plotted as functions of z at three different times 
in the Compton unit: t = l(black), t = 500 (blue) and t = 1500 (red). As shown 
in Fig. [TJ the electric field E(z, t) oscillation at the center {z = 0) is most rapid, 
and gets slower and slower at spatial points {\z\ > 0) further away from the center. 
This implies the electric field wave propagating in the space, and the directions of 
propagations are indicated. 

region TZ, the radiation fields (I25|26l) approximately are 



N 
LU 



-0.5 




5 



E rad (£,x) 



c 2 |x 



1 



/ 



dt' 



- ret 



(27) 
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Fig. 3. The net charge density t) [see Eq. ([15]) ] as a function of 2 at three 
different times: t = 1 (black, nearly zero), t = 500 (blue) and t = 1500 (red). It is 
shown that the net charged density value \p(z, t)\ is zero at the center where the 
initial electric field gradient vanishes [see Eq. (I20j) ]. whereas it increases as the initial 
electric field gradient increases for \z\ > 0. 
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Fig. 4. The net electric charge density p(z, t) [see Eq. (j 1 5 1) ] as a function of t at three 
different points: z = (red, nearly zero) , z = ij\ (blue) and z = 1/2 (black). It is 
shown that the net electric charge density p(z, t) (except the center z = 0) increases 
as time. 



where we use the charge conservation (|9]) and total neutrality condition of 
pairs J n d 3 x.'p(t', x') = 0. The first terms in Eqs. ( I25f26j) are the Coulomb- 
type fields decaying away as 0(1/ |x| 2 ). The Fourier transforms of Eqs. ( 1271) 
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Fig. 5. Electric current densities j z (z,t) [see Eq. (|14p ] as functions of z at three 
different times: t = 1 (black), t = 500 (blue) and t = 1500 (red). Following Eq. (fT9|) . 
the electric current alternates following the alternating electric field (see Fig. [1]), 
the plateaus indicate the current saturation for v ~ c and its spatial distribution is 
determined by the initial electric field E ext (z). 
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Fig. 6. Electric current densities j z (z,t) [see Eq. (fbij) ] as functions of t at three 
different points: z = 1/2 (black), z = £/4 (blue) and z = (red). The plateaus 
(see also Fig. [6]) for the current saturation values increases as time, mainly due to 
the number-densities n± of electron-positron pairs increase with time. In addition, 
they are maximal at the center z = where the initial electric field is maximal, and 
decrease as the initial electric field E ex t(z) decreasing for \z\ > 0. 

and fl2HD are 



E rad (cu,x) 



c 2 |x| 



-D(w), B rad (w, x) w R x E rad (w, x) 



D(w) = / d 3 x / dt'e% 
Jn Jt 



<9J(t',x') 
dt' 



(29) 
(30) 



where the wave number k = u/c and the numerical integration (!30|) is carried 
out overall the space-time evolution of the electric current J(x', £') (see Figs. [6] 
and El). For defmiteness we thinks of the oscillation currents occurring for some 
finite interval of time T or at least falling off for remote past and future times, 
so that the total energy radiated is finite, thus the energy radiated per unit 
solid angle per frequency interval is given by 



d 2 I 

2|D(u;)| 2 . (31) 



dudVt 



The squared amplitude |D(w)| 2 as a function of uj gives the spectrum of the 
radiation (see Fig. EJ), which is very narrow as expected with a peak locating 
at Wpeak ~ 0.08m e = 4KeV for Eq = E c , consistently with the plasma oscil- 
lation frequency (see Fig. [TJ. The energy-spectrum and its peak are shifted 
to high-energies as the initial electric field-strength increases, and the relation 
between the spectrum peak location and the electric field-strength is shown 
in Fig. [HJ In addition, the energy-spectrum and its peak are also shifted to 
high-energies as the temporary duration T of plasma oscillations increases 
(see Fig. [TJ. In calculations, the temporary duration T = 3500tc< is chosen, 
not only to satisfy the adiabaticity condition Eq. (123]) [39] . but also to be in 
the time duration when the oscillatory behavior is distinctive (see Figs. [HUH]), 
since the oscillations of pair-induced currents damp and pairs annihilate into 
photons [22J. The radiation intensity ( 13T1) depends on the strength, spatial 
dimension and temporal duration of strong external electric fields, created by 
either experimental setups or astrophysical conditions. 
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Fig. 7. In the Compton unit, normalizing D{uj) [see Eq. ()3Q|) ] by the volume 
V = f <i 3 x' of the radiation source J(i',x'), we plot \D{oj)\ 2 [see Eq. (EIJ)] rep- 
resenting the narrow energy-spectrum of the radiation field E ra d and peak locates 
at the frequency w pC ak ~ 0.08m e . 
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Fig. 8. The peak frequency w pe ak of the radiation approximately varies from 4KeV 
to 70 KeV as the initial electric field strength Eq varies from E c to 10-E C . The 
values for very large field-strengths Eq/E c > l possibly receive corrections, since the 
semiclassical pair-production rate (j22j) is approximately adopted and the pressure 
term (see [32]) is not properly taken into account. 



Conclusions and remarks. We show the space and time evolutions of pair- 
induced electric charges, currents and fields in strong external electric fields 
bounded within a spatial region. These results imply the wave propagation of 
the pair-induced electric field and wave-transportation of the electromagnetic 
energy in the strong field region. Analogously to the electromagnetic radia- 
tion emitted from an alternating electric current, the space and time varia- 
tions of pair-induced electric currents and charges emit an electromagnetic 
radiation. We show that this radiation has a the peculiar energy-spectrum 
(see Fig. [7]) that is clearly distinguishable from the energy-spectra of the 
bremsstrahlung radiation, electron-positron annihilation and other possible 
background events. This possibly provides a distinctive way to detect the ra- 
diative signatures for the production and oscillation of electron-positron pairs 
in ultra-strong electric fields that can be realized in either ground laboratories 
or astrophysical environments. 

As mentioned in introduction, the critical electric field E c will be reached 
soon in ground laboratories and sensible methods to detect signatures of pair- 
productions become important. Recently, the momentum signatures of pair- 
production is found (10] in a time- varying electric field E(t) with sub-cycle 
structure. On the other hand, space-based telescopes the Swift-BAT |JT], NuS- 
TAR jl2] and Astro-H [33] focusing high-energy X-ray missions, will also give 
possibilities of detecting X-ray radiation signature, discussed in this paper, 
from compact stars with electromagnetic structure. 
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